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Both  solid  to  gas  (CaS04  +  SiO  — T CaSiOj  +  S02  and  Si3N4  +  3C  (di¬ 
amond)  — *  and  3SiC  -f  2  Nj)  can  be  shock-induced  at  ballistic  velocities. 
Because  of  the  endothermic  nature  of  the  gas-producing  reactions,  the  extent 
of  reactions  observed  are  limited  to  interfaces.  We  found  the  above  reactions 
proceed  to  a  much  less  extent  than  calculation  by  equilibrium  thermodynamic 
calculations.  Reaction  products  are  found  to  be  10-J  times  those  calculated  for 
equilibrium.  We  show  that  the  extent  of  reaction,  rather  than  limited  by  the 
usual  diffusion  processes  appears  to  be  controlled  by  dynamic  mixing  processes 
arising  from  Rayleigh-Taylor  instabilities  at,  for  example,  CaSi04-Si02  inter¬ 
faces.  We  apply  a  theory  developed  by  Drucker  [1]  to  account  for  the  observed 
extent  of  reaction. 
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SEM  micrograph  of  shot  1100.  Light  area  is  S^fy,  voids  where  5 
diamond  grains  plucked  out  indicate  weak  bonding. 


2:  SEM  micrograph  of  shot  1106.  Darker  area  is  quartz,  lighter  5 

area  anhydrite. 

3a:  SEM  micrograph  of  shot  1099  near  sample  edge.  Dark  gray 

area  is  SiO,  light  gray  area  is  "Ca4Si3SOi3.  Black  area  is  epoxy.  5 

3b:  Shot  1099  away  from  edge.  Large  grains  are  SiO,  fine-grained  5 

material  is  CaSCX}. 

4.  X-ray  diffraction  spectra  of  silica-anhydrite  mixture  before- 
and  after-shock  (Shot  #1109,  42.3  GPa).  The  two  spectra  at 
bottom  are  JCPDS  standards  for  CaS04  and  Si02-  Molar 
ratio  of  silica  to  anhydrite  is  3:1.  Initially  crystalline  quartz 

is  amorphized  in  the  after-shock  material.  9 

5.  Atomic  ratio  of  Ca/O  profiles  of  three  recovery  shots  on 
calcium  sulfate  and  quartz.  A  cold  pressed  sample  was  used 
to  estimate  SEM  resolution.  Theoretical  fits  are  given  by 
Equation  12  with  parameters:  (1)  dashed  curve:  1  =  0.53  pm, 

L  =  0;  (2)  solid  curve:  1  =  0.53  pm,  L  =  1.5  pm.  The  data  at 

1  pm  on  the  dashed  curve  represent  three  analyses  at  different 
locations  near  the  boundary.  1 1 

6.  Atomic  ratio  of  S/Ca  profiles  of  recovery  shots  on  calcium 
sulfate  and  quartz.  Solid  curve  is  given  by  Equation  14. 

Further  into  SiC>2,  S/Ca  ratio  becomes  indeterminate  and  is 
indicated  by  the  dashed  line.  Error  bars  represent  SEM 
analytical  uncertainty.  Results  for  other  20  mm  shots  are 

similar  and  not  shown.  12 

7.  Degree  of  devolatilization  via  reactions  (A)  and  (B)  as  a 
function  of  shock  pressure.  The  curves  are  calculated  using 

Equations  16-21.  Partial  pressures  of  O2  and  SO2/SO3  are 
taken  to  be  0.2  and  10"4  bar,  respectively,  in  the  calculation, 
which  are  representative  values  for  normal  atmosphere. 

Impact  velocities  are  indicated  for  shock  pressures  of  20, 

40,  60,  and  80  GPa.  19 


lafck 


Page 

3 


1:  Recovery  experiments  on  three  gas-producing  reactions 

2:  Equation-of-state  constants  of  anhydrite,  silicate  and 

their  mixture.  7 
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Abstract 


Both  solid  to  gas  (CaSO.,  +  SiO  — *  CaSi03  4-  SO>  and  SisN^  +  3C  (di¬ 
amond)  — *  and  3SiC  +  2  N2)  can  be  shock-induced  at  ballistic  velocities. 
Because  of  the  endothermic  nature  of  the  gas-producing  reactions,  the  extent 
of  reactions  observed  are  limited  to  interfaces.  We  found  the  above  reactions 
proceed  to  a  much  less  extent  than  calculation  by  equilibrium  thermodynamic 
calculations.  Reaction  products  are  found  to  be  10~'’  times  those  calculated  for 
equilibrium.  We  show  that  the  extent  of  reaction,  rather  than  limited  by  the 
usual  diffusion  processes  appears  to  be  controlled  by  dynamic  mixing  processes 
arising  from  Rayleigh-Taylor  instabilities  at.  for  example.  CaSiO^-SiOj  inter¬ 
faces.  We  apply  a  theory  developed  by  Drucker  [1]  to  account  for  the  observed 
extent  of  reaction. 
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1  Introduction 


Shock-induced  endothermic  solid-solid  phase  changes  and  gas-producing 
reactions  have  been  studied  in  material  consolidation  and  synthesis,  espe¬ 
cially  starting  with  powdered  materials  (e.g.[2.  3.  4.  5]).  Although  ther¬ 
mal  effects  induced  by  shock  compression  is  the  major  factor  inducing 
chemical  reactions  under  dynamic  stress  loading,  unloading  within  the 
time  scale  of  microseconds  can  be  expected  to  proceed  very  differently 
than  under  equilibrium  conditions.  Gas-producing  reactions,  as  well  as 
solid-solid  phase  changes  involving  large  density  changes  with  their  poten¬ 
tial  effects  on  partitioning  of  linear  momentum,  may  have  some  ballistic 
application.  Particularly,  we  have  examined  the  following  three  reactions: 

CaSO-i(anhydrite)  -f  SiOv  — ♦  CaSi();d wollastonite)  +  SO,*  j  (A) 
CaSO,  +  SiO  —  CaSiQa  +  SO,  T  (B) 

SiaXT  -I-  3C  (diamond)  — *  3SiC  +  2 XX  f  (C) 

2  Gas-producing  reactions 

A  series  of  recovery  experiments  and  bulk  chemical  analyses  were  con¬ 
ducted  on  the  recovered  samples.  Xo  (solid)  reaction  products  are  found 
for  reactions  (1)  and  (2).  but  reaction  (3)  is  discovered  to  take  place  to 
a  varied  extent  in  different  shots.  These  results  contradict  Gibbs’  forma¬ 
tion  energy  calculations,  which  allow  all  three  reactions  to  proceed  below 
calculated  after-shock  temperatures. 

The  shock  experiments  utilize  a  20  mm  gun  at  Caltech's  Shock  Wave 
Lab.  Reactants  (in  powder  form,  particle  sizes  10-30  //in.  except  Si^.Xb  is 
whisker  shaped)  are  mixed  and  pressed  in  304  stainless  steel  containers. 
The  sample  chambers  are  usually  evacuated  to  30  millitorr  until  just 
before  the  shots  are  fired.  The  shot  parameters  are  listed  in  Table  1.  Flyer 
plates  are  tantalum  unless  specified  otherwise  in  parentheses.  Hngoniot 
pressures  are  calculated  following  the  formulation  outlined  by  Yang  *  t  <//. 
[5]  assuming  lOO'/f  crystal  density,  actual  reflected  pressures  are  higher, 
but  lion-uniform. 
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Table  1:  Recovery  experiments  on  three  gas-producing  reactions 


Using  reactant  Hugoniots  as  release  paths  (assuming  the  main  part  of 
the  reaction  takes  place  during  the  late  stage  of  release),  the  calculated 
average  after-shock  temperatures  are  2700-3000°C:  for  shot  1100.  because 
of  the  inverted  mass  ratio,  the  same  temperature  for  shot  1106  is  1200T. 
Recovered  samples  are  analyzed  with  SEM.  electron  microprobe,  and  x- 
rav  diffraction.  The  SEM  micrographs  of  shots  for  which  reaction  (1) 
was  studied  (Fig.  1)  show  SUN.}  is  molten  after  shock,  but  diamond  is 
intact.  Although  in  some  small  regions  near  boundaries  of  the  two  reac¬ 
tants  decrease  in  nitrogen  I\u  x-ray  intensity  relative  to  silicon  intensity 
is  detected,  it  is  believed  to  be  an  effect  of  geometric  absorption  by  neigh¬ 
boring  carbon  atoms.  All  prominent  peaks  in  XRD  spectra  are  identified 
with  dianiond/SiaX^  peaks,  therefore  the  reaction  product  must  be  lower 
than  the  detection  limit  (about  2c/c ).  and  gas  production  is  not  enough  to 
account  for  the  sample  holder  explosions  reported  by  Yang  ef  al.  For  shot 
1066.  reaction  (2)  does  not  occur,  but  an  interesting  feature  is  discovered 
in  the  sample:  quartz  grains  apparently  became  molten  but  anhydrite 
did  not  (Fig.  2)  despite  SiO/s  higher  melting  point  ( 1 500‘C  vs.  CaSO/s 
HOOT). 

In  comparison  to  the  inert  reactions  (1)  and  (2).  reaction  (3)  is  quite 
active.  The  bulk  reaction  product  (CaSiOO  yield  is  estimated  at  about 
301/  for  shot  1098.  Reduction  of  a  factor  of  4-5  in  sulfur  content  is 
seen  while  calcium  is  unchanged,  and  the  deficit  in  atomic  number  is 
approximately  made  up  with  incorporation  of  silicon  which  demonstrates 
the  importance  of  SiO  participation.  To  find  the  initiation  of  the  reaction, 
shot  1099  was  conducted  at  a  lower  pressure.  Clearly  separated  reaction 
and  no-reaction  zones  are  seen  in  the  SEM  micrographs  (Figs.  3a.  b). 
Both  materials  appear  to  have  been  molten.  In  the  no-reaction  zones, 
there  may  be  some  melting  in  SiO.  especially  at  grain  boundaries,  but 
CaSO^  remains  solid  during  compression.  Like  quartz.  SiO  also  has  a 
higher  melting  point  (  1700°C)  than  CaSO^. 


2.1  Reaction  between  silica  and  anhydrite 


Fig.  1:  SEM  micrograph  of  shot  1100. 
Light  area  is  S&N4,  voids  where  diamond 
grains  plucked  out  indicate  weak 
bonding.  TJA93133SFD 


Fig.  3a:  SEM  micrograph  of  shot  1099 
near  sample  edge.  Dark  gray  area  is  SiO, 
light  gray  area  is  'CaiSiaSOn.  Black 
area  is  epoxy.  TJA93135SFD 


Fig.  2:  SEM  micrograph  of  shot  1106. 
Darker  area  is  quartz,  lighter  area 
anhydrite. 

TJA93134SFD 


Fig.  3b:  Shot  1099  away  from  edge. 
Large  grains  are  SiO,  fine-grained 
material  is  CaSQ|. 

TJA93136SFD 
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Experimental  parameters  are  listed  in  Table  1  (two  shock  experiments 
on  ankvdrite/SiO  are  listed  and  will  be  discussed  later).  The  starting 
material  was  a  mixture  of  silica  (crystalline  or  amorphous.  Alfa  #13024 
and  89709)  and  natural  crystalline  anhydrite  (Ward’s  Geology  #46E0535) 
powders  and  was  pressed  into  target  container  to  initial  densities  of  60- 
85%  of  its  Archimedian  density.  Average  silica  grain  size  was  4  pm. 
and  anhydrite  grains  were  mostly  between  30-100  pm.  Equation-of-state 
constants  for  the  mixture  (see  Table  2)  were  calculated  from  pre\ious 
anhydrite  data  of  Simakov  et  a/. [8]  and  quartz  data  of  Swegle  ct  a/.[9] 
using  the  formulae  by  Boslough  [2]:  assuming  uniform  stress  distribution, 
for  a  two-component  system. 


Ks 


V"  =  miV‘j  4- 

A  05  =  [(l*l/A*0Sl)  +  (*‘2 /A  052)]" 

A'o2sh(l  +  K'0Sl )/ A’q5j  +  r2(l  +  K'QS2)/K‘U 


-  1. 


(1) 

(2) 

(3) 


where  /«,,  1'.  r I\0Si  and  K'os,  are  the  components'  mass  fractions,  spe¬ 
cific  volumes,  initial  volume  fractions,  bulk  moduli  and  their  pressure 
derivatives  at  zero  pressure.  Shock  pressures  determined  by  impedance- 
match  method  range  from  27.4  to  42.3  GPa.  Five  20  mm  shots  and  the 
two  40  mm  shots  (see  Table  1)  all  had  1:1  molar  ratio  of  anhydrite  : 
silica.  The  initial  porosity  of  the  mixture  were  17.2%  (110C)  to  40.0% 
(1108):  Shot  1107  employed  fused  quartz:  The  40  nun  shots  were  con¬ 
ducted  to  determine  the  shock  duration  effect  on  the  reaction,  but  no 
reaction  was  seen  in  shot  917.  and  shot  923  was  not  recovered:  Finally, 
shots  1109  and  1110  had  anhydrite  :  silica  molar  ratios  of  1:3.7  and 
3.O.T.  Recovered  samples  were  analyzed  with  petrographic  microscopy, 
scanning  electron  microscopy  (SEM,  instrument:  Camscan  Series  2  with 
Tracor  Northern  EDS  detector  TH-3/54-6901,  operated  at  15  kY)  and  X- 
ray  diffraction  (XRD.  instrument:  Scintag  DMC-008.  radiation  source: 
Cu-I\Ul).  Compared  with  the  original  material,  the  changes  exhibited 
in  the  20  mm  (except  shot  1111.  which  will  be  described  separately  in 
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Table  2.  Equation-of-state  constants  of 


anhydrite,  silica  and  their  mixture. 


Material 

P  (g/'-i*n 

A'os  (GPa) 

Aos 

Anhydrite  (LPP) 

2.97 

38.5 

6.0 

Silica  (HPP) 

4.29 

350 

3.3 

Mixture  (1:1  molar) 

3.28 

48.7 

7.6 

the  following  paragraph)  post-shock  samples  are  quite  similar:  in  agree¬ 
ment  with  previous  research  [10],  silica  becomes  amorphous  in  spite  of 
its  original  crystallinity  (see  XRD  spectra  in  Figure  4:  Anhydrite  is  re¬ 
covered  as  a  crystalline  phase.  Although  shock-induced  mosaicism  in  the 
crystal  grains  was  observed  with  cross-polarized  light  on  petrographic 
microscope,  it  appears  unlikely  that  anhydrite  recrystallized  from  a  melt 
as  no  rounding  of  the  grains  was  observed  (see  Figure  2). 

In  shot  1111  (in  which  10  mass%  iron  powder  was  intentionally  mixed 
in  addition  to  anhydrite/quartz),  devolatilization  was  much  more  exten¬ 
sive  than  the  rest  of  the  shots  and  reaction  of  iron  to  iron  sulfate  and 
iron  sulfide  were  observed.  In  their  study  of  sulfur  speciation  in  basaltic 
glasses  [11],  Carroll  and  Rutherford  reported  that  proportion  of  dissolved 
sulfur  present  as  sulfate  (as  opposed  to  sulfide)  increases  from  near  0V£  at 
FMQ  (fayalite-magnetite-quartz)  oxygen  fugacity  to  near  100c  at  2  to 
3  log /o,  units  above  FMQ.  The  oxygen  fugacity  present  in  these  recov¬ 
ery  experiments  was  well  into  the  sulfate  stable  regime.  The  presence  of 
iron  sulfide  led  us  to  believe  the  greater  degree  devolatilization  of  CaSO^ 
within  1  mm  of  the  stainless  steel  container  in  all  the  20  mm  shots  was 
affected  bv  the  reducing  effect  of  the  metal  and  would  not  have  occurred 
in  its  absence.  In  the  central  metal-free  region,  the  dimensions  of  possible 
reaction  zones  are  so  limited  that  they  were  nearly  at  the  limit  of  spatial 
resolution  of  the  SEM.  In  the  following  three  sections  we  will  attempt  to 
derive  the  actual  degree  of  devolatilization  from  experimentally  observed 
chemical  compositions  at  different  locations  in  the  samples. 


2.1.1  SEM  instrument  resolution 


The  SEM  electron  beam  spot  is  much  less  than  1/mi,  but  the  dimen¬ 
sion  of  excitation  volume  in  the  sample,  and  therefore  the  instrument 
resolution,  is  larger  due  to  electron  scattering  and  secondary  fluorescence 
in  the  sample  [12].  A  “smearing"  function  is  assumed  to  convolve  with 


8 


25 


Quartz 


Diffraction  angla  26 


Fig.  4.  X-ray  diffraction  spectra  of  silica-anhydrite  mixture  before-  and  after¬ 
shock  (Shot  #1109,  42.3  GPa).  The  two  spectra  at  bottom  are  JCPDS 
standards  for  CaSO^  and  Si02-  Molar  ratio  of  silica  to  anhydrite  is  3:1. 
Initially  crystalline  quartz  is  amorphized  in  the  after-shock  material. 
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the  “true"  chemical  composition  to  give  the  observed  composition.  The 
function  form  is  taken  to  be  Gaussian: 

/(*)  =  ~^jex  p(-z2//2).  (4) 

r  f{x)dx  -  1  (5) 

J  —  X 

where  2/  is  a  measure  of  spatial  resolution.  Anhydrite  and  quartz  disks 
of  1  mm  thickness  each  were  sandwiched  together  and  heated  at  573  I\ 
for  6  hours,  followed  by  a  24-hour  press  at  ^  4000  psi  so  that  plastic  flow 
may  take  place  to  produce  a  good  contact  (with  less-than-1  pm  gap)  as 
a  no-reaction  reference.  A  comparison  between  the  SEM  analysis  across 
the  pressed  boundary  and  calculation  (convolution  of  Equation  4  and  a 
step  function)  found  1  —  0.53pm  to  provide  the  best  fit  (Figure  5).  This 
agrees  with  our  expectation  that  the  resolution  distance  is  larger  than 
the  electron  beam  diameter. 

The  data  in  Figures  5  and  6  are  corrected  for  secondary  fluorescence 
excited  by  characteristic  radiations.  Another  concern  has  been  the  flu¬ 
orescence  excited  by  the  continuous  spectrum.  We  use  an  approximate 
equation  ((15.10)  in  Reed  [13]).  modified  for  compounds  by  mutiply- 
ing  the  ratio  of  mass  attenuation  coefficients  of  the  excited  element  .4 
(p^-)  and  the  compound  (pc),  the  intensity  of  fluorenscence  //  relative  to 
electron-excited  characteristic  K-radiation  of  element  .4  in  the  compound 
1$  is  thus 

4  =  9.7  x  i0-*Z't£.  (C) 

k-  « 

(Z  is  the  atomic  number  of  the  excited  element).  The  correction  factor 
for  continuum  fluorescence  is 

f/  =  i/u+4>-  <■> 

Jc 

We  consider  three  cases: 
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Fig.  5.  Atomic  ratio  of  Ca/O  profiles  of  three  recovery  shots  on  calcium  sulfate 
and  quartz.  A  cold  pressed  sample  was  used  to  estimate  SEM  resolution. 
Theoretical  fits  are  given  by  Equation  12  with  parameters:  (1)  dashed  curve:  1 
-  0.53  pm,  L  -  0;  (2)  solid  curve:  1  -  0.53  pm,  L  -  1.5  pm.  The  data  at  1  pm 
on  the  dashed  curve  represent  three  analyses  at  different  locations  near  the 
boundary. 
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Fig.  6.  Atomic  ratio  of  S/Ca  profiles  of  recovery  shots  on  calcium  sulfate  and 
quartz.  Solid  curve  is  given  by  Equation  14.  Further  into  SiOj,  S/Ca  ratio 
becomes  indeterminate  and  is  indicated  by  the  dashed  line.  Error  bars 
represent  SEM  analytical  uncertainty.  Results  for  other  20  mm  shots  are 
similar  and  not  shown. 


1.  Compound  is  50  mol.%  CaSC>4  and  50  mol.%  Si02  (in  the  middle  of 
the  mixing  zone),  lfa/l£a= 1.21%.  If/Ic= 0.24%; 

2.  Compound  is  near  100%  Si(>2  with  a  trace  of  Ca  and  S  (deep  into 
silica).  lfa/lSa= 0.67%,  7///£=0.13%: 

3.  A  sharp  CaS04-Si02  boundary,  with  the  electron  beam  shifted  to 
SiC>2  side  so  that  the  electron-excitation  volume  is  completely  in  SiC>2. 
Any  Ca  and  S  signal  is  purely  due  to  secondary  fluorescence  from 
the  S1O2  continuum  (Si  characteristic  line  is  not  energetic  enough 
to  excite  Ca  or  S).  The  secondary  fluorescence  intensity,  relative  to 
electron-excited  radiation  in  pure  CaS04.  is  given  by: 

T -  =  0.5  x  9.7  x  10-8Z4^c--£&.  (8) 

IcoSCU  HCaSO< 

Equation  8  is  very  similar  to  Equation  6  except  the  factor  0.5.  which 
arises  because  only  half  of  the  continuous  radiation  goes  into  CaSCAj. 
neglecting  the  finite  width  of  the  primary  X-ray  source.  The  ratios 
for  Ca  and  S  calculated  are  0.23%  and  0.11%. 

More  detailed  numerical  calculations  were  done  and  the  results  agree 
within  ±0.3%.  Case  (3)  agrees  very  well  with  observations  at  1  pm  in 
quartz  from  the  cold-pressed  boundary,  where  Ca  and  S  signal  intensities 
are  0.2±0.1%  and  0.4±1%  of  those  in  pure  CaS04  (Figure  5).  For  the 
shock-recovered  samples,  the  Ca  and  S  intensities  within  ~  3pm  from  the 
boundaries  are  much  higher  than  the  secondary  fluorescence  level  and  the 
corrections  are  negaligible  compared  to  the  analytical  uncertainty. 

2.1.2.  Mixing  of  sulfate  with  silica 

In  shocked  samples  the  boundary  layer  between  calcium  sulfate  and 
silica  is  thicker  than  the  cold-pressed  edge.  In  the  following  we  examine 
several  possible  mixing  mechanisms: 
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1.  Solid  state  diffusion: 

The  diffusion  constants  of  H,  1|50  and  3uSi  atoms  in  quartz  have  been 
documented  in  [14].  At  800°C,  they  vary  over  a  wide  range,  with  H 
having  the  highest  D  =  2.5  x  10” 11  m2/s.  and  30 Si  having  the  lowest 
D  =  1.3  x  10-21  m2/s.  In  the  time  scale  of  our  experiments  ~  1/is. 
the  characteristic  distance  y/Df  ~  10_8-10“2/im  is  much  smaller  than 
the  observed  reaction  zone  thickness: 

2.  Liquid  state  diffusion: 

Rubie  el  al.  directly  measured  oxygen  self-diffusivitv  in  N^Si^Oy 
melt  up  to  1825°C  and  between  4-10  GPa  [15].  The  diffusion  con¬ 
stant  they  reported  ranges  from  1.0  to  4.2  x  10" 10  m2/s  increasing  with 
temperature  and  pressure.  Si-0  bond  breaking  is  the  basic  process 
controlling  both  O  self-diffusion  in  N^S^Oy  and  CaS04  diffusion  in 
silica  melt  which  is  of  present  interest.  It  is  possible  to  use  these 
data  to  obtain  an  order-of-magnitude  estimate  of  mixing  time  and 
length  scales  in  the  Si02  liquid.  Again,  in  the  1/is  shock  duration, 
the  highest  diffusivity  (4.20xl0_1°  m2/s)  gives  characteristic  distance 
2  x  10-2//m,  which  is  still  too  small  to  account  for  the  pm-size  mixing 
layer; 

3.  Rayleigh-Taylor  instability: 

Although  initially  crystalline  quartz  is  amorphized  during  shock,  we 
cannot  conclude  that  it  has  been  once  molten  because  quartz  can 
transform  to  diaplectic  glass  without  melting  [10].  Since  there  is  a 
strong  contrast  in  strength  of  quartz  and  anhydrite  (~  1  GPa  for 
quartz  and  ~  0.1  GPa  for  anhydrite),  we  suggest  Rayleigh-Taylor 
instability  as  the  third  mixing  mechanism. 

Rayleigh-Taylor  instability  arises  at  interfaces  between  two  materials 
of  different  strength  when  they  are  strongly  accelerated  or  decelerated 
along  a  direction  perpendicular  to  their  planar  interface.  According 
to  the  theoretical  model  by  Drucker[l].  when  shock  wave  propagates 
from  the  stronger  material  into  the  weaker  material,  the  interface 


14 


experiences  alternating  compressional  and  tensile  stress  due  to  per¬ 
turbations  (bumps)  on  material  surface.  When  the  stress  difference 
exceeds  the  strength  of  the  stronger  material  (<7o).  the  bumps  grow 
freely  and  instability  occurs.  Two  important  derivations  of  the  theory 
are  threshold  perturbation  amplitude: 

C  =  //(I  +  jt/2  )a„F/P.  (9) 

and  time  dependence  of  instability  growth  above  the  threshold: 

h  -  htf  =  (h0  -  h$)cost\j3Pl/(\pH).  (10) 

where  H  is  material  thickness,  ho  is  initial  perturbation  amplitude.  /> 
is  material  density.  P  is  shock  pressure.  A  is  the  perturbation  wave¬ 
length.  F  is  a  geometric  factor  between  0.25-1.  and  3  =  ^—7  ~  2-  S. 

A  prominent  feature  of  the  theory  is  that  the  threshold  is  independent 
of  wavelength  A.  Experiments  by  Barnes  et  al.  on  aluminum  and  304 
stainless  steel  plates  support  the  theoretical  prediction  [16). 

For  our  case,  <70  ~  1  GPa,  p=2.65  g/cm3.  P=40  GPa.  A  H  grain 
size  4/im.  the  threshold  thickness  perturbation  calculated  from  Equa¬ 
tion  9  is  ~  0.4/mi.  which  is  very  reasonable.  Above  threshold,  the 
growth  is  very  fast:  for  an  initial  perturbation  10%  above  threshold. 
Equation  10  indicates  that  it  takes  a  few  nanoseconds  to  grow  to  sizes 
comparable  to  A/2,  after  which  the  theory  is  no  longer  valid. 

From  the  above  discussion,  Rayleigh-Taylor  instability  emerges  as  a 
plausible  mechanism  to  drastically  increase  surface  area  at  contact  in¬ 
terface  and  at  the  same  time  reduce  the  distance  material  has  to  diffuse 
to  mix  efficiently.  We  thus  propose  the  reaction  to  be  Rayleigh-Taylor 
instability-driven,  and  only  takes  place  at  a  later  stage  at  solid-solid  in¬ 
terfaces.  Such  reactions  would  have  the  following  characteristics: 

1.  The  amount  of  reaction  is  surface- limited  and  therefore  independent 
of  grain  size,  molar  ratio,  etc.  of  the  reactants.  This  is  consistent 
with  our  observation: 
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2.  Rayleigh- Taylor  instability  as  a  mechanical  process  does  not  rely  on 
compositional  gradient,  therefore  reaction  rates  have  a  less  important 
role  in  these  reactions  compared  to  in  those  driven  by  diffusion. 


Although  we  propose  the  anhydrite- Si  O2  reaction  is  not  dominantly  a  dif¬ 
fusive  process,  we  mathematically  describe  the  shocked  sample  interface 
by  a  linear  diffusion  profile.  Denote  gca(x)  =the  molar  ratio  of  Ca/O.  the 
solution  of  one  dimensional  diffusion  equation  (heat  equation,  see.  e.g.. 
[17])  with  initial  conditions: 


0Co(*)|t=O  = 


0.25 

0 


if  x  <  0 
if  x  >  0 


is 

9ca(*)  =  ^ertc(r/L).  (11) 

A  fit  to  the  shocked  sample  profile  (convolved  with  f{x))  gives  the 
mixing  length  L  ~  1.5/im  (solid  line  in  Figure  5): 

9ca W ( * )  -  /( -r  ~  *' )9Ca (x')dx'  (12) 


2.1.3  Degree  of  devolatilization 


The  experimentally  seen  devolatilization  shows  some  scatter  (Figures  5. 
6).  For  simplicity,  we  assume  an  average  of  20%  S  molar  loss  in  the  region 
-I  <  jt  <  L.  so  the  ratio  S/O  is 


^(j)  =  0.8pca(*r)  =  0.1erfc(j/I). 


The  S/Ca  profile  observed  on  SEM  is  given  by 


(S/Ca)(  j*) 


/*x  gcaW)?~{r‘-r)i/,'dr' 


(13) 


(14) 
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Bulk  devolatilization,  defined  as  the  fraction  of  sulfur  loss  in  the  bound¬ 
ary  layer  over  the  total  sulfur  mass  in  the  original  sample,  is  given  by 

DY  =  -  x/-t(gCa  -  (15) 

9Ca  I* - XR 

where  R  is  the  anhydrite  grain  size.  The  factor  3  takes  into  account  the 
three  dimensional  effect. For  R  —  100pm,  evaluation  of  the  formula  yields 
a  numerical  value  of  DV  =  6  x  10“3. 

We  infer  that  anhydrite  mixed  into  silica  upon  shock  loading  anti  un¬ 
derwent  devolatilization  during  release.  Tvburczy  and  Ahrens  [7]  used  the 
following  approach  tc  calculate  the  extent  of  shock-induced  reactions: 

1.  Entropy  excess  required  for  incipient  reaction: 

Sm  =  £"  'fdT.  (1C) 

where  Tq  is  room  temperature.  Tir  is  the  temperature  of  incipient 
reaction  (at  which  the  sums  of  Gibbs  formation  energies  for  the  re¬ 
actants  and  products  are  equal),  Cp  is  the  atmospheric  pressure  heat 
capacity  at  constant  pressure. 

2.  Entropy  excess  required  for  complete  reaction: 

Scr  =  Sir  +  AS  -  Y,  H.Phi(P/Po).  (17) 

|U  frojwri* 

where  AS  is  the  entropy  difference  between  reactants  and  products 
as  computed  from  Robie  et  ai  [6],  the  last  term  on  the  right  takes 
into  account  effects  of  partial  pressures,  P,.  of  gas  products  (P0  is  the 
ambient  pressure),  n,  is  the  number  of  moles  of  gas  specie  i.  R  is  the 
gas  constant  (8.31  J  mol_1K_1). 

3.  Entropy  gain  in  the  shocked  state  (and  in  the  post-shock  state  as¬ 
suming  isentropic  release): 

±SH  =  Str  +  CMTH/Ts)'  (IS) 

where  Str  is  the  entropy  change  of  phase  transition  during  compro*- 
sion.  T$  is  the  temperature  of  isentropic  compression  from  initial 
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volume  Vo  at  temperature  T,  to  Hugoniot  volume  Vh  and  isentropic 
pressure  Ps : 


T,  =  7*  exp 


(19) 


where  *y  is  the  Griineisen  parameter.  Th  is  the  shock  temperature, 
determined  by: 


~(Ph  -Ps)=  C“  C,dT 

*>  ■'T, 


(20) 


C,  is  heat  capacity  at  constant  volume  (at  high  pressure). 


4.  The  extent  of  reaction  is  given  by: 


Fraction  of  material  reacted  = 


-  Sjr 

Scr  -  Sift 


(21) 


The  results  from  Equation  21  for  reaction  (A)  are  shown  as  the  solid 
and  broken  curves  in  Figure  7.  The  equation-of-state  parameters  are  the 
same  as  in  Table  2.  and  5<r=0  for  the  calculation.  Local  devolatilization 
in  the  reaction  layer  is  close  to  theoretical  calculation,  but  decomposed 
anhydrite  is  only  a  small  portion  (~  6x  10~3)  of  the  total  mass  because  of 
the  dimension  of  the  reaction  layer  is  much  smaller  relative  to  the  grain 
size  (~  100//m). 

To  examine  the  role  of  enthalpy  difference  in  shock-induced  reactions, 
we  conducted  two  recovery  experiments  on  anhydrite  and  amorphous 
silicon  monoxide  (SiO,  Alfa  #89430)  powder  mixture  (Table  1).  The 
reaction  in  question  is 

CaSO*  +  SiO  —  CaSi03  +  SO?  T  (B) 

We  observed  extensive  reaction  in  the  recovered  material  of  shot  1098. 
where  S  in  anhydrite  is  reduced  by  a  factor  of  4-5.  Although  enthalpy  of 
SiO  glass  is  unknown,  it  is  a  less  stable  compound  than  quartz  and  there¬ 
fore  absorbs  less  energy  in  reaction  (B).  The  veaker  bonding  structure 
compared  to  SiO*  could  also  give  rise  to  a  mucn  higher  diffusion  rate. 
Bo\h  may  contribute  to  the  excessive  reaction  wlmh  occurred. 
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Fig.  7.  Degree  of  devolatilization  as  a  function  of  shock  pressure  of 
reaction  CaSCXi+SiOs  — *  CaSi03+S(>3  The  curves  are  calculated 
using  Equations  16-21.  Partial  pressures  of  0;>  and  SO3  are  taken  to  be 
0.2  and  10”4  bar  respectively  in  the  calculation,  which  are  representative 
values  for  normal  atmosphere.  Impact  velocities  are  indicated  for  shock 
pressures  of  20.  40.  60  and  80  GPa. 
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3  Conclusions 

We  studied  three  gas-producing  solid-solid  reactions  under  shock  condi¬ 
tions: 

CaS04 (anhydrite)  +  Si02  — ♦  CaSiOsfwollastonite)  +  SO3  T  (A) 
CaS04  +  SiO  — ♦  CaSiOa  4-  SO2  T  (B) 

Si3N4  +  3C  (diamond)  —  3SiC  +  2X2  T  (C) 

The  reactions  are  found  to  proceed  only  in  a  very  limited  reaction 
layer  at  the  interfaces  of  the  corresponding  reactants,  and  for  this  reason, 
reacted  portion  of  the  material  is  a  very  small  percent  (~  10“2)  of  the 
bulk,  although  inside  the  reaction  zone,  the  extent  of  reaction  can  be 
much  greater  (~  209c)  and  close  to  theoretical  calculations. 

With  emphasis  on  reaction  (A),  we  propose  these  reactions  are  con¬ 
trolled  by  the  dynamic  process  of  Rayleigh-Taylor  instability  instead  of 
diffusion,  as  in  reactions  under  equilibrium  conditions.  A  model  devel¬ 
oped  bv  Drucker  is  applied  to  anhydrite-SiO-2  interface.  Since  diffusion 
and  reaction  rates  are  secondary  restraints  for  these  reactions,  such  a 
mechanism  explains  the  lower  degree  of  reaction  compared  to  equilib¬ 
rium  calculations  observed  in  all  three  reactions. 
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